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In  this  report,  we  present  results  of  an  analysis  of  the  isoplanatic 
aspects  of  predetection  compensation  imagery.  We  have  developed 
formal  theoretical  results  for  how  the  achievable  modulation 
transfer  function  for  predetection  compensation  varies  with  the 
angular  size  of  the  image,  or  with  the  angular  separation  between 
the  image  and  the  reference  source.  Detailed  evaluation  is 
carried  out  for  propagation  over  a  path  with  uniform  optical 
strength  of  turbulence  along  the  path  (nominally  a  horizontal 
path).  It  is  found  t'  *it  for  propagation  along  a  horizontal  path 
of  length  z,  if  the  wavelength  and  strength  of  turbulence  are  such 
as  to  result  in  a  coherence  length  rQ,  then  the  angular  separation 
between  the  image  and  the  reference  must  be  less  than  5Q  *  r0/2z 
if  predetection  compensation  is  to  produce  useful  imagery  of  fine 
details . 

In  order  to  extend  the  calculations  of  the  isoplanatic  aspects  of 
predetection  compensation  to  include  results  for  ground-based  \ 

viewing  of  space  objects,  it  is  necessary  to  have  a  model  for  the 
vertical  distribution  of  the  optical  strength  of  turbulence.  We 
have  generated  such  a  model,  making  use  of  the  low  altitude  air¬ 
craft  thermal  probe  measurements  of  Koprov  and  Tsvang,  and  the 
I  high  altitude  balloon  thermal  probe  data  of  Bufton.  The  model  is 
presented  in  terms  of  a  table  of  numerical  values,  along  with  a 
discussion  of  the  data  processing  details  by  which  the  model  was 
extracted  from  the  available  raw  data. 

As  part  of  the  problem  of  measuring  the  optical  effects  of  turbu¬ 
lence  over  finite  non-horizontal  propagation  paths,  such  as 
between  an  aircraft  and  the  ground,  this  report  considers  the 
possibility  of  using  differential  angle  of  arrival  measurements. 

The  advantage  of  this  type  of  measurement  is  that  it  allows  tracking 
of  a  point  source  and  mrasures  wavefront  distortion  in  terms  of 
(differential)  angle  of  arrival  at  two  sub-apertures  on  the  same 
mount,  giving  meaningful  results  without  placing  significant  con¬ 
straint  on  the  tracking  accuracy  of  the  mount.  Significantly,  the 
basic  equipment  for  this  type  of  measurement  is  available  at  RADC. 
The  questions  to  be  answered  by  theoretical  analysis  are  how  does 
the  mean  square  differential  angle  of  arrival  vary  with  separation 
of  the  two  subapertures,  and  is  the  expected  value  for  the  largest 
achievable  separation  great  enough  compared  to  equipment  measurement 
precision  to  allow  for  meaningful  measurements.  These  questions 
are  answered  in  this  report.  The  dependence  is  evaluated  and 
tables  of  results  are  presented.  It  is  found  that  measurement 
accuracies  are  adequate  to  allow  propagation  path  characterization 
with  modest  (i.e.,  5%  to  1050  precision. 
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Summary 


This  report  presents  the  more  significant  and  complete  results 
to  date  on  a  program  of  theoretical  study  related  to  unconventional  imagery 
and  techniques  for  suppression  of  atmospheric  turbulence  effects  on  image 
quality.  The  objective  of  this  work  has  been  to  develop  an  understanding 
and  quantitative  theory  of  those  aspects  of  propagation  through  atmospheric 
turbulence  that  would  limit  the  performance  of  unconventional  imagery 
techniques,  and  to  develop  a  basis  for  quantifying  the  effective  magnitude 
of  atmospheric  turbulence. 

In  this  report,  we  present  results  of  an  analysis  of  the  isoplanatic 
aspects  of  predetection  compensation  imagery.  We  have  developed  formal 
theoretical  results  for  how  the  achievable  modulation  transfer  function  for 
predetection  compensation  varies  with  the  angular  size  of  the  image,  or 
with  the  angular  separation  between  the  image  and  the  reference  source. 
Detailed  evaluation  is  carried  out  for  propagation  over  a  path  with  uniform 
optical  strength  of  turbulence  along  the  path  (nominally  a  horizontal  path). 

It  is  found  that  for  propagation  along  a  horizontal  path  of  length  z  ,  if  the 
wavelength  and  strength  of  turbulence  are  such  as  to  result  in  a  coherence 
length  tq  ,  then  the  angular  separation  between  the  image  and  the  refer¬ 
ence  must  be  less  than  i?0  =  r0/2z  if  predetection  compensation  is  to 
produce  useful  imagery  of  fine  details. 

In  order  to  extend  the  calculations  of  the  isoplanatic  aspects  of 
predetection  compensation  to  include  results  for  ground-based  viewing  of 
space  objects,  it  is  necessary  to  have  a  model  for  the  vertical  distribution 
of  the  optical  strength  of  turbulence.  We  have  generated  such  a  model, 
making  use  of  the  low  altitude  aircraft  thermal  probe  measurements  of 
Koprov  and  Tsvang,  and  the  high  altitude  balloon  thermal  probe  data  of 
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Bufton.  The  model  is  presented  in  terms  of  a  table  of  numerical  values, 
along  with  a  discussion  of  the  data  processing  details  by  which  the  model 
was  extracted  from  the  available  raw  data. 

As  part  of  the  problem  of  measuring  the  optical  effects  of  turbu¬ 
lence  over  finite  non-horizontal  propagation  paths,  such  as  between  an 
aircraft  and  the  ground,  this  report  considers  the  possibility  of  using 
differential  angle  of  arrival  m  asurements.  The  advantage  of  this  type 
of  measurement  is  that  it  allows  tracking  of  a  po’nt  source  and  measures 
wavefront  distortion  in  terms  of  (differential)  angle  of  arrival  at  two  sub¬ 
aperture  j  on  the  same  mount,  giving  meaningful  results  without  placing 
significant  constraint  on  the  tracking  accuracy  of  the  mount.  Significantly, 
the  basic  equipment  for  this  type  of  measurement  is  available  at  RADC. 

The  questions  to  be  answered  by  theoretical  analysis  are  how  does  the 
mean  square  differential  angle  of  arrival  vary  with  separation  of  the  two 
subapertures,  and  is  the  expected  value  for  the  largest  achievable  sepa¬ 
ration  great  enough  compared  to  equipment  measurement  precision  to 
allow  for  meaningful  measurements.  These  questions  are  answered  in 
this  report.  The  dependence  is  evaluated  and  tables  of  results  are  pre¬ 
sented.  It  is  found  that  measurement  accuracies  are  adequate  to  allow 
propagation  path  characterization  with  modest  (  i.  e.  ,  5%  to  10%)  precision. 
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ABSTRACT 


In  Chapter  I,  the  concept  of  predetection  compensation  imagery  is 
defined  and  a  formal  analysis  of  the  effects  of  atmospheric  turbulence  on 
the  compensated  image  MTF  is  set  up.  An  expression  is  obtained  for  the 
MTF  in  terms  of  an  integral  over  the  propagation  path.  The  integrand  is 
a  function  of  the  ref  r a ctive  -  index  structure  constant,  the  image  frequency, 
f  ,  the  angular  separation  between  the  target  object  and  the  reference 
source,  i?  ,  and  the  angle,  0  ,  between  the  vector  aspects  of  f  and  0 

For  horizontal  propagat’on,  it  is  found  that  two  parameters  can  be  defined 
which  govern  the  quality  of  the  compensated  MTF.  These  are  d0  ,  the 
predetection  compensation  critical  angle,  and  fQ  ,  the  transition  frequency 
We  find  that  for  image  frequencies  below  f0  ,  there  is  no  useful  compen¬ 
sation,  in  fact  we  suffer  an  effective  doubling  of  the  strength  of  the  turbu¬ 
lence.  (However,  if  [Q  is  small  enough,  this  is  of  no  consequence.)  For 
image  frequencies  greater  than  fQ  ,  the  MTF  is  down  from  diffraction- 
limited  by  no  more  than  a  factor  of  exp  [-2.  58  (,?/i?0f/3]  .  So  long  as  & 


is  made  less  than  i?„  ,  the  compensation  will  be  quite  effective.  Expres- 


dn  and  f_  are 


*o  = 


rQ/2z  and  f  =  2i?z/\  ,  where 


sions  obtained  for  „0  ami  i0-«,z/A  ,  wnere  z 

is  the  path  length  and  rQ  is  the  critical  coherence  length  for  conventional 
long  exposure  imagery. 


In  Chapter  II,  medium  and  high  altitude  thermal  probe  data  by 
Bufton,  and  low  altitude  thermal  probe  data  by  Koprov  and  Tsvang  are 
processed  to  provide  a  smooth  distribution  for  the  refr  a  ctive  -  index  struc¬ 
ture  constant  in  the  atmosphere.  The  model  is  used  to  calculate  the 
critical  length  for  diffraction-limited  imaging,  i.  e.  ,  r0  ,  and  the  log- 
amplitude  var iance ,  i.  e.  ,  a^2  .  Comparison  with  measured  values  shows 
very  good  agreement  for  r0  ,  and  only  fair  agreement  for 


al 
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In  Chapter  III,  the  problem  of  measuring  differential  angle  of  arrival 
through  two  relatively  small  apertures  of  variable  separation  is  considered. 
The  motivation  is  to  evaluate  the  practicality  of  a  measurement  program 
using  such  a  quantity  to  observe  atmospheric  turbulence  wavefront  distor¬ 
tion  effects  when  one  end  of  the  link  is  moving  and  can  not  be  tracked  pi 
cisely.  Theoretical  results  for  the  mean  square  difference  in  angle  of 
arrival  are  developed.  Numerical  results  have  been  calculated  using  a 
computer.  Comparison  of  expected  magnitude  of  effects  with  available 
measurement  instrument  precision  indicates  that  the  experiment  should 
be  possible,  but  will  depend  on  our  ability  to  achieve  an  rms  single  axis 
angle  of  arrival  measurement  precision  of  the  order  of  0.  1  arc  seconds 
or  better. 
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CHAPTER  I 

Isoplanatic  Aspects 
of 

Predetection  Compensation  Imagery 


1 


Introduction 


The  concept  of  isoplanatism  refers  to  the  fact  that  over  some  size 
field-of-view  the  image  of  an  object  will  suffer  the  same  degradation  from 
diffraction-limited  quality,  independent  of  where  the  object  is  within  that 
field-of-view.  It  is  common  usage  to  refer  to  the  field-of-view  as  the 
isoplanatic  patch.  As  thus  defined,  the  concept  of  isoplanatism  refers 
basically  to  the  performance  of  a  lens  or  mirror  imaging  system.  How¬ 
ever,  this  same  concept  can  be  applied  in  discussing  the  effect  of  atmos¬ 
pheric  turbulence  on  an  otherwise  ideal  imaging  system.  It  is  in  this 
context  that  we  shall  consider  the  subject  of  isoplanatism. 

Before  proceeding  farther,  we  wish  to  make  certain  disclaimers 
and  place  certain  bounds  on  the  extent  of  the  problem  we  shall  be  treating. 
First  of  all,  we  call  attention  to  the  fact  that  it  is  by  no  means  clear,  or 
even  likely  that  atmospheric  turbulence  induced  resolution  effects  will  be 
strictly  constant  over  any  non-trivial  size  field-of-view.  Rather,  it  is 
likely  that  whatever  resolution  variation  there  is  over  the  field-of-view 
will  show  a  monotonic  dependence  on  the  size  of  the  field  angle.  We  could 
define  the  isoplanatic  patch  as  being  determined  by  the  size  of  the  field 
angle  that  makes  the  variation  just  reach  some  (arbitrarily  chosen)  level 
of  significance.  However,  rather  than  engage  in  a  matter  of  such  arbitrary 
nature,  we  shall  simply  seek  knowledge  of  the  field-angle  dependence, 
which  we  shall  refer  to  as  the  isoplanatic  dependence  of  resolution,  and 
avoid,  as  much  as  possible,  discussion  of  the  size  of  the  isoplanatic  patch. 

Secondly,  we  wish  to  take  note  of  the  fact  that  there  are  a  variety 
of  distinct  conventional  and  unconventional  imaging  techniques  presently 
being  considered  for  imaging  through  atmospheric  turbulence.  Each  of  these 
techniques  functions  in  a  different  manner  and  is  sensitive  to  somewhat  dif¬ 
ferent  aspects  of  the  optical  effects  of  atmospheric  turbulence.  There  is, 
therefore,  no  reason  to  believe  that  the  isoplanatic  dependence  of  any  two 
of  these  techniques  will  be  the  same.  In  this  paper,  we  shall  be  considering 


one  of  the  more  interesting  of  the  unconventional  imaging  techniques, 
which  we  refer  to  as  predetection  compensation.  There  is  no  assur¬ 
ance  that  the  isoplanatism  dependence  that  we  shall  develop  in  this  paper 
for  predetection  compensation  imaging  will  be  applicable  to  any  other 
imaging  technique,  and  we  caution  the  reader  against  assuming  that  this 
result  necessarily  is  relevant  to  any  of  these  other  techniques.  (At  a 
later  time,  we  hope  to  extend  our  work  to  cover  some  of  these  other 
techniques.  ) 

Predetection  Compensation 

Predetection  compensation  is  the  generic  term  used  to  describe 
any  one  of  a  number  of  techniques  that  seek  to  estimate,  in  real  time, 
the  nature  of  the  atmospheric  turbulence  induced  wavefront  distortion 
at  the  imaging  system's  collection  aperture,  and  in  real  time  compensate 
for  the  wavefront  distortion.  In  this  way,  the  image  as  formed  is  com¬ 
pensated  for  turbulence  effects.  (The  term  "predetection  compensation" 
is  perhaps  best  understood  in  juxtaposition  to  the  term  "postdetection 
compensation.  "  Postdetection  compensation  is  the  process  of  correcting 
the  distorted  image  after  it  is  formed,  by  some  form  of  modulation  trans¬ 
fer  function  compensation,  or  by  some  more  or  less  equivalent  technique.  ) 

The  key  to  predetection  compensation  lies  in  two  areas.  The  first 
of  these  is  the  ability  to  in  real  time  estimate  the  nature  of  the  phase  dis¬ 
tortion  and  perhaps  also  the  intensity  variation  at  the  aperture.  The  second 
key  is  the  ability,  in  real  time,  to  modify  the  optical  system  so  as  to  cor¬ 
rect  for  the  phase,  and  possibly  also  the  intensity  variations,  so  that  the 
corrected  optical  signa’l  contains  only  the  information  descriptive  of  the 
object  being  viewed.  Techniques  are  under  development  by  several  organi¬ 
zations  aimed  at  demonstrating  both  of  these  capabilities  and  producing  a 
predetection  compensated  image. 

It  would  be  beside  the  point  to  go  into  the  details  of  the  techniques 
that  are  being  considered.  Our  concern  here  is  with  the  evaluation  of  one 


potential  limitation  in  the  predetection  compensation  concept.  This  limita¬ 
tion  relates  to  the  isoplanatic  aspect  of  the  process.  We  are  concerned 
with  the  fact  that  in  order  to  estimate  the  phase  and  intensity  variations, 
all  techniques  that  have  been  proposed  require  the  observation  of  a  refer¬ 
ence  source.  The  estimated  variation  is  then  used  to  correct  the  wave 
coming  from  some  point  on  the  object  of  interest.  In  general,  the  field 
angle  associated  with  the  reference  source  will  not  be  identical  to  the 
field  angle  for  the  point  on  the  object  we  are  imaging.  Our  concern  is  to 
develop  a  quantitative  understanding  of  how  the  angular  separation  between 
the  reference  source  and  the  point  being  imaged  affects  the  resolution  of 
the  compensated  image.  We  consider  this  to  be  the  problem  of  isoplanatisrr 
for  predetection  compensated  imagery. 

Analytic  Problem  Definition 

For  the  purposes  of  analysis,  we  shall  consider  the  problem  of 
predetection  compensated  imaging  of  a  point  target  source  at  field  angle 
"^l  ,  using  a  point  reference  source  at  field  angle  ~Q3  .  We  shall  study 

the  ensemble  average  modulation  transfer  function  for  the  predetection 
compensated  target  image,  ( Troc  (?))  ,  at  image  frequency  T  (in  cycles 

rer  radian  field-of-view),  as  a  function  of  the  angular  separation 
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In  many  propagation  situations,  intensity  variations  are  minor, 
i.  e.  ,  the  log-amplitude  variation  is  much  less  than  one  neper,  and  so 
there  is  no  great  utility  in  correcting  for  intensity  variations.  Only  phase 
variation  correction  would  be  of  concern.  In  other  propagation  situations, 
however,  intensity  variations  could  be  sizable,  and  it  would  be  necessary 
to  utilize  intensity  compensation  in  order  to  obtain  a  well  compensated 
image.  While  such  compensation  is  somewhat  more  difficult  than  phase 
compensation,  there  is  nothing  fundamental  that  prevents  such  compen¬ 
sation  from  being  utilized.  In  all  cases,  the  use  of  intensity  compensation 


along  with  phase  compensation  will  give  at  least  as  good,  if  not  a  signi¬ 
ficantly  better  compensated  image  than  would  the  use  of  phase  compen¬ 
sation  alone.  Consequently,  to  investigate  the  ^-dependence  of  the  pre¬ 
detection  compensated  image,  we  shall  treat  the  case  in  which  both  phase 
and  intensity  variations  are  corrected. 

We  shall  use  the  notition  and  0(x.'Q1)  to  denote  the  instan¬ 

taneous  turbulence  induced  variation  of  log  -  amplitude  and  phase,  respec¬ 
tively,  at  a  point  x  on  the  collection  aperture  plane,  for  a  wave  coming 
from  the  point  target  source  at  ^  .  Similarly,  we  shall  use  ^(x,^)  and 
0  (x,  03)  to  denote  the  same  quantities  for  the  wave  coming  from  the  point 
reference  source  at  Qs  .  In  predetection  compensation,  we  form  estimates 
°f  9S)  and  0(x,l^a)  ,  which  we  shall  assume  for  the  purpose  of  this 

analysis  are  perfect  estimates,  and  correct  the  received  wave  from  the 
target  accordingly.  This  means  that  we  form  our  image  using  a  wave  with 
residual  log-amplitude  and  phase  errors  at  x  ,  given  by 

A*(x,£)  =  jUx.'Bj)  -  l(x,~§3)  ,  (2a) 

A0(x,£)  =  0(x,'01)  -  0(x,~93)  ,  (2b) 

respectively.  If  there  were  no  field  angle  dependence,  so  that  ^(x,^) 
exactly  equaled  £(x,  9a)  and  0  (x,  0 x)  exactly  equaled  0(x,^a)  ,  then  the 
residual  errors  would  be  zero  and  the  predetection  compensated  image 
would  be  Ideal.  However,  to  the  extent  that  isoplanatism  is  less  than 
perfect,  the  residual  errors  will  be  non-zero  and  the  predetection  com¬ 
pensated  image  will  be  less  than  ideal. 

The  calculation  of  the  modulation  transfer  function  now  proceeds 
in  almost  exactly  the  same  manner  we  have  utilized  previously1  for  calcu¬ 
lation  of  the  long  exposure  MTF  for  conventional  imagery.  If  we  let  U(x) 
denote  the  scalar  form  of  the  compensated  electromagnetic  field  at  aperture 
position  x  ,  and  let  u(y)  denote  the  associated  image  scalar  form  of  the 
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electromagnetic  field  at  the  position  y  in  the  focal  plane,  then  it  is  well 
known  that 


u(y)  =  A  f  dx  U(x)  exp  (-ikx  •  y) 


(3) 


where  A  is  a  constant  of  proportionality  and  k  =  2tt/x  is  the  optical  wave 
number.  (Here  and  throughout  this  paper  integrals  without  explicitly  stated 
limits  are  to  be  understood  as  being  definite  integrals  with  infinite  limits. 

In  Eq.  (3),  the  effective  limits  of  the  integration  are  provided  by  the  fa<’t 
that  U(x)  vanishes  for  x  outside  the  limits  of  the  collection  aperture.  ) 
The  instantaneous  power  distribution  in  the  focal  plane  is  proportional  to 
|  u(y ) |  2  .  Since  the  MTF  is  just  the  fourier  transform  of  the  impulse 

response,  which  |u(x)|2  represents,  we  can  write  for  the  MTF  at  image 
frequency  T 


tpdc  (?)  =  B  r  dy  J  u(y)[2  exp  (2tt  i  T  •  y)  ,  (4) 

where  B  is  a  normalization  constant  chosen  to  make  Tpp^O)  =  1  . 

If  we  combine  Eq.  's  (3)  and  (4),  we  get 

tpoc^)  =  A2B  JJJ  dy  dx  dx'  U*(x')  U(3c)  exp  |~2niy  •  (?+  ^  .  (5a) 


Noting  that  carrying  out  the  y-integration  gives  rise  to  a  delta  function 
which  then  allows  the  x'-integr ation  trivially,  we  get 


Troc(f)  =  A2B  Jdx  U*(x  -  X?)  U(x) 


(5b) 


Now  introducing  W(x)  to  define  a  circular  aperture  of  diameter  D  , 
according  to  the  equation 


W(x)  =  / 


1  , 


if 


*  *  D 


0  ,  if  x  >  £  D 


(6) 
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we  can  write  the  compensated  scalar  field  at  the  aperture,  U(x)  in 
terms  of  the  compensated  log -amplitude  and  phase  errors  bl(x,d)  and 
A0  £,£),  respectively.  We  have 

U(x)  =  W(x)  exp  [A£(x,j?)  +  i  A{>  (x,i?)]  ,  (7) 

and  if  we  substitute  this  into  Eq.  (5b),  we  get 

t „(.(?)  =  A2B  J  dx  W(x  -  \?)  W(x)  exp  {[A£(x,i?)  +  Ajfc(x-X7,jj>)] 

+  i  [ A0 (3c)  -  A0(x-\T,#)]}  .  (8) 

The  ensemble  average  predetection  compensated  modulation  transfer 
function  can  now  be  written  as 

<  TP0C(?))  =  A2  B  J  dx  W (x-\T)  W (x)  (exp  {[A£(x,i?)  +  A£(x-\?, ,?)] 

+  i  [A0(x,tf)  -  A0(x-\TJ)]}>  ,  (9a) 

Since  the  wavefront  distortion  statistics  are  homogeneous  so  that  the 
ensemble  average  on  the  right  hand  side  of  Eq.  (9a)  is  not  actually  a  func¬ 
tion  of  x  ,  we  can  write 

<Tpoc^)>  =  x0L  (f)  M(XTJ)  ,  (9b) 

where  TDL(f)  is  given  by  the  expression 

tol (f)  =  A2  B  J*  dx  W(3c-\?)  W(3c)  ,  (10) 

and  is,  as  we  shall  see,  the  MTF  of  a  diffraction-limited  circular  aperture 
of  diameter  D  .  M(\f*,i?)  is  given  by  the  expression 


M(XT.?)  =  (  exp  {[A£(x,»?)  +  AX  (x- \?,  i? )]  +  i  [A0(x,t?)  -  A0(x-\i\  j?  ]  }  )  .  ()1) 


The  integral  on  the  right  hand  side  of  Eq.  (10)  can  be  recognized  as  the  area 
of  overlap  of  two  circles  of  diameter  D  whose  centers  are  displaced  a 
distance  Xf  .  Making  use  of  a  little  trigonometry,  it  is  easy  to  show  that 

-  f  {  cos-1  ft)  -  ft)  [,  -  ft)J]  ,  in 


where  we  have  anticipated  the  fact  that  M(0,$)  is  equal  to  unity  and  adjusted 
the  normalization  constant  B  so  that  the  MTF  will  have  unity  value  at  f  =  0  . 
We  recognize  rDL(f)  in  Eq.  (12)  as  being  the  well  known  expression  for  the 
MTF  of  a  diffraction-limited  optical  system  with  a  circular  aperture  of 
diameter  D  . 

Our  problem  at  this  point  is  reduced  to  evaluation  of  M(\lf,i?)  . 

To  the  extent  that  this  quantity  equals  unity,  the  predetection  compensated 
modulation  transfer  function  is  ideal,  and  to  the  extent  that  it  deviates  from 
unity  (we  shall  see  that  it  is  always  less  than  or  equal  to  unity),  the  pre¬ 
detection  compensated  modulation  transfer  function  is  less  than  ideal.  The 
balance  of  the  analysis  in  this  paper  is  devoted  to  the  evaluation  of  M(X?,$)  . 
In  the  next  section,  we  develop  an  expression  for  M  in  terms  of  quantities 
describing  optical  propagation  statistics,  and  in  the  section  after  that,  we 
treat  the  propagation  problem  to  obtain  a  formulation  of  these  statistics. 

Representation  for  M(X?,|?) 

In  order  to  reduce  the  right  hand  side  of  Eq.  (11),  we  start  by 
noting  that  if  a  and  0  are  independent  gaussian  random  variables,  each 
with  zero  mean  value,  and  if  a  and  b  are  arbitrary  constants,  then  it  is 
easy  to  show  that 
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(exp  (aa  +  b0)>  =  exp  [£  (a2(a2)  +  b2  <32))] 


(13) 


Because  the  log-amplitude  i  and  the  pha.'c  q)  are  themselves  gaussian, 
it  follows  that  the  compensated  log  -  amplitude  error  Li  and  the  compensated 
phase  error  A0  are  each  gaussian.  It  follows  from  consideration  of  Eq.  's 
(2a  and  b)  and  the  fact  that  the  propagation  statistics  are  stationary  in  terms 
of  the  field  angle  dependence  (in  this  case,  and  3S  )  that  Li  and  A0 
both  have  zero  mean  value.  It  can  be  shown,  using  the  fact  that  the  propa¬ 
gation  statistics  are  isotropic,  that  the  quantities  [A^fx.j?)  +  A4(x-X7,$)] 
and  [  A0(x,  j?)  -  A0(x-x7,j9)]  are  independent.  Obviously,  then,  Eq.  (13) 
can  be  applied  to  the  evaluation  of,  the  right  hand  side  of  Eq.  (11).  We  get 

M(\f,^)  =  exp  ([A4(x,tf  )  +  A£(7-X?,£)]s  ) 

-  £  <[A0(x,£)  -  A0(x-x7,i?)]2  >)  .  (14) 

In  order  to  reduce  the  terms  in  the  exponent,  we  introduce  the  two 
structure  function  related  statistical  quantities 

0*#)  -  <[X(x,  Qj)  -  A(x%  91)][£(x,  02)  -  i(x,qs)-])  ,  (15a) 

and 

®0(p,?)  =  <L*a.V  -  0(x',91)][0(7,-e3)  -  0(^,f3)]>  ,  (i5b) 


where 


P  =  x  -  x' 


(16) 


We  shall  refer  to  these  two  quantities  as  the  log- amplitude  and  the  phase 
hyperstructure  functions.  (In  the  limiting  case  of  f  =  it,  so  that  £  =  0 
the  hyperstructur e  functions  reduce  to  the  ordinary  structure  functions.  ) 
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If  we  make  use  of  Eq.  (2b),  we  see  that  the  phase  dependence 
in  Eq.  (14)  can  be  written  as 

<[A0(x,i?)  -  A0(x  -  \?,i?)]a) 

=  <  ([0(x,'e1)  -  0(x,e2)]  -  [0(x-\?, ex)  -  0(x-\f\"e2)]f  > 

=  <{[0(x,e1)  -  0(x-x?,?1)]  -  [0(x,  e2)  -  0(^-x?,‘4)])a  > 

=  2  [$0(X?,O)  -^(\f,£)]  .  (17) 

Reduction  of  the  log-amplitude  dependence  in  Eq.  (14)  is  somewhat  more 
complex.  It  is  most  conveniently  accomplished  by  introducing  the  log- 
amplitude  hyper  covariance  function 

Sx(p,3)  =  ([Afx.fy  -  I][£(^,0a)  -  7])  ,  (18) 

where 

1  =  U(x,e)>  .  (19) 

It  is  easy  to  show  that 

<  f[i(x,?1)  -  i(3Ua)]  ±  [  l(x'  -  £(x',  l^)]}3  > 

=  4  {[®A(0,  0)  -  £x(0,tf)]  ±  [Ejfp.  0)  -  £  EA(p,?) 

-*$*(?.-?>]}  •  (20) 

By  making  use  of  Eq.  (20),  we  can  write  the  log-amplitude  dependence  of 
the  exponent  in  Eq.  (14)  as 


<[Af(x,i?)  +  Af(x-X?,i?)]3  ) 

=  (ftiGc.fy  -  A(x,es)]  +  [^(x-xT.eJ  -  l (x-X?."3j>)]}2  > 


(21) 


=  -  <{[£(x,"e1)  -  a(x,?2) j  -  [x(x-x^,"ej)  -  x(x-x?, 9a)]}s ) 

+  8  [S£(0,0)  -  Ex(0,^)] 

=  -  <  r  [  X  (x ,  )  -  £(x-x7',"0l)]  -  [£(x,0^)  -  £(x-xf>,'§r)]}s  ) 

+  8  [<S£(0,  0)  -  <srx(0,?)] 

=  -  2[®X(X?\  0)  -  S£(x?,?)]  +  8[E£(0,  0)  -  E£( o,?)]  . 

Now  if  we  substitute  Eq.  's  (17)  and  (21)  into  Eq.  (14),  we  gee 

M(X?,0)  =  exp  {-  [%(Xf,  0)  -  SX(X?,?)]  +  4  [£^(0,  0)  -  SX(O,0)] 

-  [!ty(X?,0)  -  S0(xr,tf)]] 

=»  exp  {-  [®x(x?,0)  -  Sx(x?,?)]  -  [5)0  (xT,  0)  -  ®0(X?,?)]]  ,  (22) 

where  we  have  dropped  the  Ex(0,0)  -  £x(0,i>)  dependence  in  the  exponent 
since  it  does  not  represent  an  image  frequency  dependence.  Its  contribu¬ 
tion  is  multiplicative  and  is  absorbed  into  the  normalization  constant  B 
in  going  from  Eq.  (9a)  to  Eq.  (9b).  (As  we  shall  see,  dropping  this 
£- dependence  insures  that  M(X7,$)  will  have  unity  value  for  7  =  0.)  Eq.  (22) 
represents  the  endpoint  of  this  section  --  a  representation  of  M(x?,j?)  in 
terms  of  the  propagation  statistics.  In  the  next  section,  we  take  up  the  prob¬ 
lem  of  calculating  these  statistics. 

Propagation  Statistics 

The  problem  of  calculating  the  propagation  statistics  associated  with 
the  hyperstructure  functions  has  been  treated  previously.  2  The  published 
results  can  be  written  in  the  form 

%(p,i?)  =  ^  k  /  dv  CN3  J*  do  [1  -  exp  (i  a  •  "p)]  a-11/3 

o 

X  [cos  (2a  •  dv)  -  cos  (v  a2/k)]  ,  (23a) 
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8.  16  k2 


®0(P>#)  4  n2  J  dv  CN2  J  d a  [  1  -  exp  (i  a  •  "p)]  a"11/3 


X  [cos  (2a  •  $v)  +  cos  (v  a2/k)] 


The  a- integration  is  understood  to  be  two-dimensional,  in  the  plane  per¬ 
pendicular  to  the  z-axis.  (  CN2  ,  of  course,  denotes  the  refractive-index 
structure  constant,  and  k  =  2tt/X  is  the  optical  wave  number,  ) 

If  \/e  substitute  these  results  into  Eq.  (22),  we  see  that  we  can 

wr  ite 

=  exp  [  -  ^(xT,?)]  ,  (24') 

where 

8  16  k2  ^ 

7/((\f,d)  =  -  z  a —  J  dv  CN2  J  da  [1  -  exp  (i  a  •  X?)]  a"u/3 

o'1 

x  [  1  -  COS  (2a  •  £v)]  .  (24) 

Our  problem  at  this  point  is  the  reduction  of  &j(xT,0)  to  numerical  results. 

We  take  this  up  in  the  following  sections. 


Preliminary  Reduction  of  (\T ,$) 

Without  having  to  make  any  assumptions  as  to  the  nature  of  the 

distribution  of  CN2  along  the  propagation  path,  we  can  carry  out  the 
— ♦ 

a- integration.  To  do  this,  we  first  rewrite  Eq.  (24)  in  the  form 


2  n2  I  dv  f  da  a"U/3  f1  -  exP  (i  a  *  X?) 
o 

-  £  exp  (2  i  a  .  *?v)  -  £  exp  (-2  i  cr  .  $v) 

+  £  exp  [  i  a  •  (XT  +  2tfv)]  +  £  exp  [  i  a  •  (\?-2^v)]}  . 


(25) 


If  we  treat  the  a- integration  as  being  in  polar  coordinates  and  perform 
the  angular  portion  of  that  integration,  then  the  exponentials  each  give 
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Strictly  speaking,  this  equation  is  not  directly  applicable  to  Eq.  (27)  since 
in  our  case,  p  =  8/3  and  v  =  0  ,  so  the  condition  Re  p  >  Re  v-1  is  not 
satisfied.  There  is  also  a  problem  of  the  divergence  of  the  integral  with  1 
in  place  of  JQ  .  In  fact,  however,  because  we  are  dealing  with  a  difference 
of  factors,  the  a-integration  in  Eq.  (27)  is  convergent.  The  divergence  is 
associated  with  the  behavior  of  the  integrand  around  a  =  0  ,  for  which 
ct  ^3  JQ  (ax)  clearly  leads  to  a  divergence.  However,  the  differences  are 
such  that  the  quantity  in  the  curly  brackets  in  Eq.  (27)  goes  as  a2  in  the 
vicinity  of  a  =  0  ,  and  so  the  actual  dependence  near  a  -  0  is  a-2^3  . 

This  does  not  have  a  divergence.  This  means  that  we  can  use  Eq.  (28)  to 
evaluate  Eq.  (27)  [with  1  replaced  by  J0  (0  a)]  ,  and  by  an  argument  based 

on  the  principle  of  analytic  continuation  take  the  result  to  apply  for  the 
parameters  in  Eq.  (27).  Thus  we  obtain  the  result  that 
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(29) 


Z 

=  “p  2-8/3  p|y|y||  ^  Jdv  Cn2  l-(\tf/3  -  (2e?vf/3 

0 

+  £  [Xf)2  +  2(20v)  (xf)  cos  0  +  (20V)2]6/6 

+  [(>.f)2  -  2(20v)  (\f)  cos  0  +  (20V)3]5/6} 

Here  0  denotes  the  angle  between  0  and  T  . 

Equation  (29)  represents  the  most  general  reduction  of  results  we 
can  carry  out  without  making  some  assumption  regarding  the  distribution 
of  the  refractive-index  structure  constant  along  the  propagation  path,  i.e.  , 
the  dependence  of  C^2  on  v  .  In  the  next  section,  we  consider  the  parti¬ 
cular  case  corresponding  to  horizontal  propagation  for  which  Q,2  is  constant, 
independent  of  v  . 


Reduction  of  ^(X?,^)  for  Horizontal  Propagation 

If  we  restrict  attention  to  the  case  of  horizontal  propagation  for  which 
CN2  is  constant,  we  can  simplify  Eq.  (29)  by  taking  C^2  outside  the 
v-integration.  We  can  then  simplify  the  parametric  dependence  of  the  inte¬ 
gral  by  making  the  variable  of  integration  correspond  to  20v/Xf  .  If  we 
make  this  transformation  and  extract  all  possible  parametric  factors,  we 
obtain  from  Eq.  (29) 


tf?(\7,0) 


2-8/3  r(-5/6) 

n  r  (n/6)  K  20 


CN3  J  dx  { 1  +  x5/3 


-  |  [1  +  2x  cos0  +  x5  ]B/6  -  &  [  1  -  2x  cos  0  +  x2  J6'6  }  .  (30) 


Here  we  have  utilized  the  quantity  fp  ,  which  we  call  the  transition  image 
frequency,  to  denote 

f0  =  20  z/x  .  (31) 
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At  this  point,  it  is  interesting  to  seek  to  introduce  the  quantity 
r0  ,  which  measures  the  effect  of  wavefront  distortion  of  conventional 
imaging.  r0  corresponds  to  the  practical  upper  limit  imposed  by  atmos¬ 
pheric  turbulence  on  useful  optics  diameter  for  diffraction-limited  imagery.1 
For  infinite  plane  wave  propagation,  r0  is  defined  by  the  equation 


6.  88 

r  -6/3 
r0 

=  2.  91 

k2 

z  CN 

2 

to  write 

8.  16 

2-3/3 

r(-5/6) 

6_ 

.  88 

(Xf)^3 

TT 

Lt 

r(n/6) 

2 

.  91 

2t>  z 

8.  16 

2-8/3 

r  (-5/6) 

6. 

88  > 

'{  \3/3 

TT 

r(n/6) 

2. 

91  \ 

^  J  ' 

(32) 


/f  n8^3  /  2i9  z  \6/3 

6-88fc)  (ir)  1(f/w)  ■ 


(33) 


where 

<o/f 

I  (f/fo  ,0)  =  J  dx  { 1  +  xP/3  -  £  [  1  +  2x  cos  0  +  x2  ]B/6 
o 

-  £  [  1  -  2x  cos  0  +  x2  ]s/6  )  .  (34) 


The  evaluation  of  i?)  reduces  at  this  point  to  the  evaluation 

of  I(f/fQ,0)  .  In  general,  this  integral  can  only  be  evaluated  by  numerical 
techniques,  which  we  have  performed.  Before  presenting  the  numerical 
results,  it  is  appropriate  to  discuss  the  asymptotic  evaluation  of  l(f/fg,0) 
for  f/fQ  »  1  and  for  f/^j  «  1  .  We  shall  denote  the  asymptotic  form  by 
IagyB(f/%  i0)  »  understanding  that  the  appropriate  asymptotic  form  is  to  be 
used  in  the  regions  f/%  <  1  and  f/%  s  1 
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For  f/^  »  1  ,  x  in  the  integrants  is  everywhere  very  small,  and 

we  can  make  the  approximation 


1  +  x5/3  -  £  (1  +  2x  cos  0  +  Xs)6/8  -  £  (1  -  2x  cos  0  +  x2)6/8 


1  +  x6/3-  £  j  1  +  Yp  (2x  cos  0  +  x2)  +  ^  ~  (2x  cos  0)2  | 


S  +  yp  (-  2x  cos  0  +  x2)  +  —  (-2x  cos  0  )2 


X?/3  -  Y  X2  +  YY  x2  COS2  0 


(35) 


Carrying  out  the  integration  is  now  trivial,  and  we  get 


“  s  &-)  [i  - 1? 


if 


(36) 


For  f/lg  «  1  ,  the  value  of  the  integrand  in  Eq.  (34)  is  dominated 

by  the  values  lor  which  x  is  very  much  greater  than  unity,  and  we  can 
make  the  approximation 


1  +  x6'3  -  £  (1  +  2x  cos  0  +  x2f/6  -  £  (1  -  2x  co»  0  +  x2)5^6 

=  1  +  x^3  -  £  x6/3[(l  +  2x_1  cos  0  +  x'2f/e  +  (1  -  2X'1  cos0+x'2)5/e] 


1  +  x^3  -  £  x6/3 
5/6 


1  +  pT  (2x_1  cos  0  +  x'2)-f  ^5//6^2',1  ~  (2x_1  cos  0  J2  j 
4  x6^1  +  YP  (-2X"1  cos  0  +  x-3)  +  (_2x_1  cos  0)2 


1  -J  x-1/3  +  —  x“1/3  COS2  0 


(37) 
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Here  again,  carrying  out  the  integration  is  trivial,  and  we  get 


1 .„.(f /<,.*)«.  (4-)  [*  -f  (4;/3<1  -3  cosa«]  • 


if  f<fo  .  (38) 


We  have  written  a  computer  program  for  the  evaluation  of 
>  and  comparison  of  this  quantity  with  1.,^  (f /^>  ,  0 )  .  The 
computer  program  is  listed  in  Table  I.  2,  and  the  numerical  results  are 
given  in  Table  I.  3.  As  can  be  seen  from  these  numerical  results,  I  a  ( f ,0) 
is  a  very  good  approximation  to  I  for  all  values  of  f/%  when  cos30  is 
near  unity,  and  is  also  a  good  approximation  when  cos20  is  significantly 
different  from  unity,  if  f  >  $,  ,  or  less  than  about  0.1  ^  . 

With  these  asymptotic  forms,  we  see  that  Eq.  (33)  can  be  re¬ 
written  as 


where 


5  /  f  \1/3 

"  4  \1^J  ^  "  a  CO£?0)j  >  if  f  <  4) 

20  /L  \1/3 

~27  ("T- )  ^  ^  cos20)J  ,  if  f  5 


(39) 


The  angle  j?0  can  be  considered  to  characterize  isoplanatism  for 
predetection  compensation  imaging  over  a  horizontal  path.  We  shall  call 
it  the  critical  angle  for  predetection  compensation  imaging. 

In  the  next  section,  we  present  a  discussion  of  the  significance  of 
our  results. 
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Discussion  of  Re sul ts 


If  we  substitute  Eq.  (39)  into  Eq.  (24'),  and  that  into  Eq.  (9b), 
we  obtain  the  basic  result  that 


<TpDC(?)>  =  T0L(f)  ( 


(  5  /f  4/3 

exp  |-6.  88^—  j  |  1  - ~  (~^j  (I-5  coss0)j  f  if  f  <  ^  , 
exp  {-2.SS  (^-)  |  l-ff  (-^-)  (1-i  cos20)jj.  if  f*^. 


Recalling  definitions  given  previously,  we  note  that  tdl  (f)  is  the  diffraction- 
limited  MTF  of  the  aperture, 


,  -  z 

h  =  X  (3D 

and 

,  -  Jo. 

0  2z  *  (40) 

0  is  the  angle  between  the  direction  to  the  reference  source  and  the 
orientation  of  the  spatial  frequency,  ?  ,  of  interest. 

In  considering  the  results  in  Eq.  (41),  we  note  first  of  all  that  for 
f  small  compared  to  the  transition  frequency,  fQ  ,  the  MTF  reduces  to 

T0L(f)  exp  [  -6.  88  (f\/r0 ]  .  If  we  had  not  used  predetection  compensation, 
the  well-known  long  exposure  result  is  of  the  form  rDL(f)  exp  [-3.44  (f\/rQp] 
Thus  we  see  that  for  low  frequencies,  the  MTF  is  actually  poorer  than  in 
conventional  imagery.  This  is  apparently  because  for  these  low  frequen¬ 
cies,  the  wavefront  distortion  compensation  we  are  using  is  inappropriate. 

By  using  it,  we  are  simply  doubling  the  mean  square  wavefront  error,  be¬ 
cause  we  are  combining  two  independent  sets  of  errors.  This  means  that 
in  practice,  we  must  have  an  fQ  small  enough  so  that  (fQ  X/rQ)  is  less 
than  unity  so  that  the  low  frequencies  are  in  a  range  where  there  is  no 
significant  wavefront  distortion,  and  doubling  the  atmospheric  turbulence 
effect  will  be  of  little  consequence. 


In  order  to  get  fj,  small  enough,  we  have  to  have  an  appropriately 
small  value  for  $  .  A  bit  of  algebraic  manipulation  quickly  shows  that 

the  requirement  that  we  make  the  transition  frequency,  ,  small  enough 
is  equivalent  to  the  requirement  that  we  use  a  reference  source  at  an  angu¬ 
lar  distance  §  which  is  less  than  the  critical  angle  for  predetection 
compensation,  . 

For  image  frequency  f  greater  than  the  transition  frequency,  , 

we  see  that  the  predetection  compensation  MTF  is  equal  to  the  diffraction- 
limited  MTF,  except  for  a  factor  of  the  order  of  exp  [-2.  58  (t?/#0)^3]  . 

Obviously,  the  key  to  achieving  a  high  predetection  compensated  MTF  at 
high  image  frequencies  is  to  use  a  reference  source  at  a  displacement 
angle  &  which  is  less  than  the  predetection  compensation  critical  angle, 

• 

Fortunately,  we  get  the  same  requirement,  i.  e.  ,  &  <  ,  for 

good  performance  at  both  high  and  low  image  frequencies.  Obviously,  the 
magnitude  of  is  of  fundamental  importance,  and  it  is  interesting  to 
see  what  it  is  for  various  conditions.  First  of  all,  we  note  from  Eq.  (32) 
that  rQ  is  proportional  to  X6/6  .  Considering  Eq.  (40),  it  is  thus  appar¬ 
ent  that  j?0  is  a  function  of  wavelength,  and  also  varies  as  \6/5  .  If  we 

combine  Eq.  's  (32)  and  (40),  we  get 

*  =  (05  26,3  k2  2  c»a ) 

=  0.0923  z -8/5  (CN2)_3'B  .  (42) 

In  Table  I.  1,  we  present  some  representative  values  of  the  predetection 
compensation  critical  angle,  ,  along  with  corresponding  values  for 

the  nominal  angular  resolution  for  uncompensated  imagery,  namely  r0  , 
where 


-  =5.42  X"1/5  z3'5  (CN2)3/6 

ro 
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(43) 


1 

Table  I.  1 

Dependence  of  t?0  and  \/r0  For  Horizontal  Propagation 

X  =  0.  5  (j,m 


C„2(m-^3) 

R  =  1 

km 

R  =  3 

km 

R  =  1 

3  km 

*>0 

(p,  r  ad) 

X/r0 

(p,rad) 

^0 

(p.rad) 

X/r0 

(tirad) 

(lirad) 

\/rQ 

(lirad) 

1  X  1  O'13 

2.  53 

98.  7 

0.  437 

190.  8 

0. 0637 

393. 

3  x  lO'14 

5.  22 

47.  9 

0.  900 

92.  6 

0.  1311 

190.  8 

1  X  10‘14 

10.  09 

24.  8 

1.  74 

47.  9 

0.  253 

98.  7 

3  x  10‘15 

20.  8 

12.  04 

3.  58 

23.  3 

0.  522 

47.  9 

1  X  Kris 

40.  2 

6.  2  3 

6.  92 

12.  04 

1.  009 

24.  8 

3  x  10~16 

82.  7 

3.  02 

14.  27 

5.  84 

2.  07 

12.  04 

1  x  10-1® 

159.9 

1.  56 

27.  6 

3.  02 

4.  02 

6.  23 

It  is  apparent  from  these  results  that  for  horizontal  propagation, 
excluding  short  paths  with  relatively  light  turbulence,  the  critical  angle, 
j?0  ,  is  less  than  the  nominal  long  exposure  resolution,  \/rQ  .  This  means 

that  in  general,  for  horizontal  propagation  paths,  predetection  compensation 
can  only  be  utilized  if  a  seL -referencing  technique  can  be  utilized.  In  gen¬ 
eral,  the  distortion  of  a  resolvable  distinct  reference  source  provides  incor¬ 
rect  predetection  compensation  information  for  imaging  of  a  target  of  interest. 

For  vertical  propagation,  entirely  independent  calculations  start¬ 
ing  at  Eq.  (29)  need  to  be  carried  out  with  an  appropriate  vertical  distri¬ 
bution  for  the  refractive-index  structure  constant,  Q2  .  It  is  likely  that 
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in  that  case  a  simple  representation  involving  an  easily  defined  pre- 
detection  compensation  critical  angle,  i?0  ,  and  a  transition  frequency, 

ig  ,  such  as  we  have  found  for  horizontal  propagation  will  not  apply.  We 
intend  to  carry  out  an  analysis  of  the  vertical  propagation  case  in  a  later 
paper.  For  the  moment,  we  remark  that  some  insight  into  the  vertical 
propagation  can  be  obtained  from  Table  I.  1. 

If  we  are  to  apply  our  horizontal  results  to  the  vertical  propagation 
case,  we  must  assume  some  effective  height  for  the  turbulent  atmosphere. 
The  corresponding  value  of  Q,2  then  follows  from  the  fact  that  X/rQ  is 
of  the  order  of  5  x  10"6  rad.  If  we  assume  an  effective  height  of  1 ,  3,  or 
10  km's,  the  corresponding  values  of  &Q  are  of  the  order  of  40  ,  14  , 
and  4  ^rad.  From  this,  it  is  apparent  that  unless  the  effective  height  is 
of  the  order  of  1  km  or  less,  only  self- refer enced  predetection  compen¬ 
sation  can  be  utilized. 
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Table  1.2 


Computer  Program  Listing  for  Evaluation  of  I(f/^,0)  and 
I«»r«  (f /^>  •  The  comPuter  language  is  BASIC. 


i^  rr:  /r,,v :  =  =■. i 

J  f \  r*.  t  ~  T  /  -  —  i.*  <i  .»•  •  • 

1 7  „  ...  ,  #  #  ...•  v  1 1  ♦  r  .  •.  ,/  •••  t  t  r  r 

or-,-’  LET 

^  1  L —  1  =  1  /  7 
9^  7  »  TT  "  9  -  n  /  7 

9]  *>  L~T  <"  = 

°1  1  LET  ^  1  =  .'>/- 

9  ]  9  ]  r,*>  99-9" /9' > 

9]  9  LE"  ~  7  =  G  /  /, 

99  9  a  r'  p 

99]  r  p"  l~  =  -  ]  9  7  0  0  *r  0  9  7^ 

nATA  "j  .  1,  .  .  'tj  .  5,  .  <,  .1  J  .  r',  .  9,  1  ,  - 

'34'1  LE"  01  =  mr  (  .  c*?.  ]  /_,  j  tJop.'AiC) 
i  l"t  r==  1-01*01/3 

O-r.  ]-g-  9_]  TQ  og 

95  1  9”  f  'JT 

9C;7  - 

°'>*  *  *  9  I  *  t  G  A?;  =  ”;c*9«;'»  DEG.” 

979  rg-'  7-]  ~g  = 

97  1  n-r,T 

77O  V7'’"  7 

09  n  LET  g=57_n/.! 

29*  LET  1  =  0 
399  LET  ''=  j  7  7 

31 *  G0.c’T  0** 

39°  LET  I  1 =  K 1 
33*  F0r  L 1 = 1  TO  7 
34*  EG"  L*=  1  TO  9 
397  LET  v  =  +D 
3^0  fjfifrro  ,-,nn 
37  0  LET  I  2  =  ••'  1 
707  LET  D 

799  9  n  9r;p  a 97 

479  LE'r  t  3=7'  1 

41*  LET  I  =  I  +  (  I 1+4*19+13)* D/ 3 

497  LET  11=13 

425  IF  ”<  1  E- 9p  r-OTf]  /jPJP. 

490  IF  ”>  1"^  GOTO  4°°’ 

479  jr  ">  ]  9970  4*0 

44*  LE'r  J=0  1  *vt  ?9*(  i_99*pt  ?1  *C9) 

45*  GOTO  47  * 


r  t  r  r 
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vT  vT 


Table  2  -  Continued 


i*'*  L~T  .j  =  -’*(  1  .''3AC9/vn|  ) 

/tnr  T'I*r  ~?;xj  1/”J  i;,IJ  I/.J 

J  7  L  0  =  ~ 

iOC  Mr-,—  70 

/,'■=;  j_''t  D~  1  ^  *  D 
40r  M"'*7  L  1 
K(7c  70m0  r'j°r 
r  r.  »>  =  J  _  —  *  r-  J  +’•*'• 

]  Oi  J^7'T'  \f  ]  -r- 
AP''1  J-  ir=r  'IQ-n 
.^On  J^r"7  7]-  />,  nc  c  )  t  o 

<=.z)7  lt-7  x=  i  +?*”*c l  +”*x 

*5''  LX?  ’Cl  =-.  c*(’:l  +  AnC  (X)  t  "  )  +  1 

T'T?TT’r-rj 

A7  ^  P7QD 

6'x  xxn 
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Table  1.3 


Computer  Printout  for  I  (f , 0 )  and  ,0)  .  The 

output  is  generated  by  the  program  listed  in  Table  2. 


j  =  r*  *  o  pr?  # 


V* 

l(f/fo,0) 

!/1.« 

P.  P  1 

i ,  pppp  +  pp 

1 . CPP--P4 

1 . cc=Z-?6 

p.  377 

n  # 

s .  p  p  p  ~ + p  i 

O , CCPP_P^ 

0  ,  c  7  p.o_",<. 

P  .  77/j 

o  %  f?  O 

7,  ?23t:  +  P  1 

c,7 c  pc- p  c 

p.  7  5  0  Z  -  "  c 

P.  777 

7  %  ”  /• 

P  .  cppr;  +  p  | 

c,  7  o  c  p  _  p  C 

C,  9  0p".(7C 

A  /*»  ^ 

?  ,  p  P  '■  «  —  +  r  J 

1  .  p/jP'?-P4 

1 . 7 /•'■n-r.’t 

0  r,  ry  -N 

i 

1  .  <.***7+  J 

1 . 44r  7 -p/; 

j  ,  <,  r>  r-  _  -  /. 

<7  P  *7 

1  .  /)  P  O  P  +  «  l 

?.  4 3  4  ~  -  n /( 

p  ,  ,/j  O  c  _  p  /. 

7,000 

9,  po 

1  .  P  =  P-+  P 1 

7.  SPAP-P/, 

3  ,  c  n  7  t  _  p  /. 

OHO 

p  .  "  9 

1  .  !  1  lE  +  'M 

4. 7/jO”-P/; 

.'j  .  7  c  pz. 

7  #  nr? 0 

P.  1  P 

1  ,  p  r  p  tt  +  p  j 

4, 3pcp-Pi 

4. 7777.7/1 

000 

p  ,  o  p 

C  .  P.PPT  +  PP 

7. 4477-P3 

7. 443F-P3 

^  #  non 

P.7P 

7,  7  7ptT  +  PP 

1 . P 1 1 r-po 

1 . p 1 7--P0 

r.ooQ 

p.  /in 

0  ,  CPPC  +  PP. 

P  ,  P  67  JT_Pp 

0,  ppp-'.po 

^ ,  -  07 

p, 

0, PPP”+pp 

3,  C777_,~p 

7.  CP  ]  ••_  7 

r-  ap  / 

•  . 

■  ' " 

1  .  4  66 -;  +  <••' 

c.  =<,o---pp 

C  t  ••.  ■  7  •  ■  _  n  r> 

■ 

in 

1  .  /'PC(r:7  +  r  p 

‘  ,  p  c  7  ■ p  p 

r'.  1  3  4  O 

•  •  r 

r*  O  /*, 

« 

1  ^  O  C  p  ’?  +  p  P 

1 . 1 P47- p 1 

1  .  1  p  7  E  -  p  1 

p .  r  ''  6 

P,  ,  0  p 

1  .  1  1  1 E+PP 

1 . 45PE-P1 

1 .43  IE- PI 

p.  907 

1  .  PP 

1  ,  Pppc+  P(3J 

1 . «44E-?1 

1 . 9  93E-P. 1 

p.97  1 

p.  ppi 

<=. pppr-p l 

4.  ? 36E-P 1 

6.771E-71 

1  .  RP? 

7  ,  P  r 

3.  3  3  3  E  -  3  1 

1 . 3407+pp 

1 . 7  4  6Z+PP 

1  .  PP? 

/l  .  ”  3 

p,  ^pp-?_  P  J 

1 . 7P 17+PP 

1  .  pppyc+pp 

1  .  ppp 

5,  pp 

0,  pppir^r*  j 

p.  664E+r'P 

7.  563E  +  3P 

1  .  ppp 

4.  nr 

1 . 4  447- P 1 

3. P43E+PP 

3 . 34c'Z  +  pp 

1  .  ppp 

n  ,  o  n 

1  . 4P*  *>?  1 

7  ,  O  C  p—  +  P  P 

7,  p  c  p,  —  +  pn 

1  .  ppp 

S.f'i? 

1  .  cqour.p  i 

4. 6647+PP 

4 . 6  4. 4  Z  +  n  P 

1  .  ppp 

0  .  PP 

1.111  TE- ^  1 

c . 3°4Z+pp 

4,  3 *5 4 E  +  p ? 

P.999 

1  P.  PP 

1  .  p  p.  ?  E  -  p  1 

6.  1 3 1 7+P3 

6 . 1 3PE+7P 

P.  09? 

op ,  pp 

c  , pppp_pp 

1  .  3337+p  1 

1 . 33SZ+7 1 

p.099 

7P,  <7  17 

3. 333E- ?? 

3. 1 34E+P 1 

?.  1  95E+P 1 

p.  009 

4  P .  p  p 

0 ,  cp  p>r_pp 

3 .  p??Z+p 1 

3. PP5Z+P 1 

P.  999 

CP,  pp 

p ,  p  p  p  •?  -  p ;  o 

3.  R63E+P 1 

3. 343E+P 1 

p.  993 

(S,P  ,  PP 

1 . 44,67-  P5 

4.  7 1 3Z+P I 

4. 7PPE+P1 

P.997 

7  P  ,  p  p 

1 . 4^37-7? 

4.  c  7  77+p  J 

c . c  5  4Z  +  p 1 

p  ,  000 

OP,  PP 

1  ,  pcptr.pp 

6. 44 17+ PI 

6. 45  3E+  p 1 

P  .  7  9  9 

p  p  ,  p  p 

1  •  1  1  1  7-  0,0 

7,  7  p  p-r+  p 1 

7. 2°6Z+p 1 

r-  ^  n 

• 

r  p  .  r  P 

1  ,  PAPT _  n r> 

7.  Pr- 1  7  +  p  1 

7.  ?pi/4-p+p  ! 

p  .  ?  9  9 
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Table  3 


Continued 


O  o  nr?" 

•  ?  • 


I 

| 


f 


h 

/f 

Ufa 

K  f/4.0) 

.0) 

1  n.m 

r 

^  1 

]  .  ft  ft  ft  ^  -  o 

1 . c  0  °  7  -  -  a 

1  . 

f;  c;  p  ~  .  r  < 

n 

nn  ^ 

ft. 

o 

C  #  7  7  o—  4.  ^  J 

ft.  97 ft"- ft'. 

rt 

c;  -  7  T?  _  9  ^ 

rt 

77  1 

n 

n  p 

?.  3 3 97+°  1 

0  p  p  p  r  .  p  c 

r\ 

n  ^  _  p  ~ 

.0  9 

npr 

r  # 

**  /< 

0  ,  c  r.  -  r*  +  r  j 

c*  0 ^ /i 1 - c c 

“ 

7 1 7---G 

O  9 

ft  -7 

c, 

r-  G 

0  ,  r  0  ft  -  *  "  1 

1  .  r 

1 . 

~  p  -  •?  _  n  /; 

O  9 

**  **,  O 

c 

'’ft 

1  ,«.fA.-r  +  c[ 
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of  the  Optical  Strength  of  Turbulence  in 
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Pr  e  sentation 


It  is  the  objective  of  this  chapter  to  present  a  plausible  model  for 
the  vertical  distribution  of  the  refractive -index  structure  constant  C^2  (h) 
as  a  function  of  altitude,  h  ,  based  on  the  most  current  available  atmos¬ 
pheric  turbulence  data.  The  actual  distribution  is,  of  course,  stochastic, 
and  very  likely  varies  to  some  extent  with  geographic  location.  However, 
our  objective  here  is  merely  to  obtain  a  reasonable  model  which  can  be 
used  to  carry  out  order -of-magnitude  type  calculations  of  the  optical 
effects  of  atmospheric  turbulence.  With  this  objective  in  mind,  we  rec¬ 
ognize  that  not  only  is  it  not  possible  for  us  to  retain  the  details  of  fine 
scale  dependence  of  q2  on  h  which  might  be  observed  at  any  instant 
of  time,  but  also  it  is  neither  appropriate  nor  significant  to  do  so.  The 
detailed  atmospheric  physics  which  can  result  in  order  -  of-magnitude 
changes  in  CN2  when  h  changes  less  than  1  km ,  at  some  high  altitude, 
is  obviously  of  potential  interest  to  the  atmospheric  physicist  who  wishes 
to  study  the  mechanism  by  which  refractive-index  variations  are  generated 
in  the  atmosphere.  However,  for  those  of  us  who  merely  wish  to  carry 
out  optical  effects  calculations,  since  the  weighting  factors  are  relatively 
slowly  and  smoothly  varying  factors  of  position  along  the  propagation  path, 
it  is  only  necessary  to  develop  an  estimate  of  the  average  value  of 
in  each  altitude  regime.  Although  actual  values  of  q2  versus  altitude, 
h  ,  may  show  spikes,  as  depicted  in  Fig.  II.  1,  it  is  entirely  permissible 
for  us  to  smooth  the  data  so  long  as  we  retain  the  correct  integrated  value 
of  CN2  in  each  altitude  regime. 

The  absence  of  direct  optical  data  from  which  the  distribution  of 
CN2  versus  h  could  be  developed  for ce s  us  to  rely  on  therma  1  probe 
measurements  of  the  temperature  structure  constant,  q2  ,  from  which 
the  refractive-index  structure  constant,  CN2  ,  can  be  calculated.  For 
this  purpose,  we  have  found  it  convenient  to  utilize  the  balloon  flight  data 
gathered  by  Bufton1  for  the  altitude  regime  between  500  m  and  20,  000  m, 


and  the  airplane  flignt  data  of  Koprov  and  Tsvang2  for  the  altitude  regime 
between  50  m  and  500  m.  (It  is,  in  fact,  the  recent  availability  of  Bufton's 
high  altitude  data  that  makes  the  development  of  this  new  turbulence  distri¬ 
bution  model  practical  and  appropriate.  ) 

The  Koprov  and  Tsvang  data,  because  it  represents  level  flight  at 
each  altitude,  gives  a  reasonable  average  value  and  so  relatively  little 
data  smoothing  is  required.  Their  data,  representing  the  average  of  seven 
day  time  flights,  and  the  average  of  five  night  time  flights,  is  shown  in 
Fig.II.  2.  The  solid  lines  show  a  fit  to  the  data  according  to  the  equations 


CN2(h)  =  6.  5  x  10-14  h-1  >, 

50  <  h  <  500  m 

(day  time) 

(la) 

C*2(h)  =  \ 

1 

|  7.  0  X  10-14  h-4*'3 

1.  5  X  10-ie 

,  50<h<100m 

,  100  <  h  <  500  m 

(night  time) 

(lb) 

where  CN2  is  in  units  of  m"3'3  and  h  is  in  units  of  m  . 

The  balloon  flight  data  of  Bufton  represents  four  night  time  flights. 
Using  the  integrated  measured  values3  of  CN2  (derived  from  C^2  )  and 

smoothing  the  data  to  suppress  sharp  altitude  dependence,  we  have  calcu¬ 
lated  the  average  value  of  CN3  in  1  km  altitude  increments  centered  at 
1  km,  2  km,  3  km,  etc.  ,  up  to  20  km,  for  each  of  the  four  flights.  The 
results  are  shown  in  Table  II.  1,  along  with  the  average  value  for  all  four 
flights.  These  values  are  depicted  in  Fig.  II.  3. 

Rather  than  slavishly  follow  the  detailed  variations  of  these  data 
points,  we  have  reasoned  that  the  sharp  variations  we  see  in  Fig.  jj>  3  are 
the  result  of  only  having  four  data  sets  to  average  over.  We  assume  that 
the  real  dependence  in  a  statistical  average  sense  is  an  altitude  regime 
dependence,  and  have  therefore  utilized  the  much  smoother  curve  shown 
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to  approximate  what  we  believe  the  data  says  is  the  actual  statistical 
average  altitude  dependence.  The  corresponding  adjusted  set  of  values 
for  C^2  is  listed  in  Table  II.  2. 


We  consider  Eq.  's  (la  and  lb)  and  the  data  in  Table  II.  2  to  repre¬ 
sent  our  model  for  the  altitude  dependence  of  the  refractive -index 
structure  constant. 

It  is  interesting  to  note  that  this  model  bears  a  qualitatively 
striking  relationship  to  that  proposed  by  Hufnagel4,  particularly  in  terms 
of  the  excess  turbulence  in  the  10  km  to  15  km  altitude  range,  although  the 
values  we  have  appear  to  be  smaller  than  those  of  Hufnagel's  model.  How¬ 
ever,  this  may  merely  be  a  consequence  of  our  decision  to  smooth  over 
altitude.  (It  is  difficult  to  be  more  quantitative  in  the  comparison,  be¬ 
cause  of  the  problem  of  estimating  the  altitude  range  for  the  disturbed 
layers  in  Hufnagel's  model.  ) 

As  a  check  on  our  model  for  C^2  ,  we  have  calculated  the  value 

of  the  critical  coherence  length,  r0  ,  and  the  log-amplitude  variance, 

<j  2  ,  for  X  =  0.  55  g,m  and  vertical  propagation,  according  to  the  equations 


and 


p*th 


(2) 


a*  =  0.56  (-^)7/S  J  dh  h^  c;3(h)  .  (3) 

P»th 

We  calculate  that  r0  should  equal  0.  095  m  in  the  day  time  and  0.  100  m 
at  night.  Calculations  of  a^2  give  essentially  the  same  result  for  both 
day  and  night,  namely  0.  116  neper2  .  The  night  time  r0  result  is  to  be 
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compared  with  measured  values6  centered  around  0.  114  m  .  This  agree¬ 
ment  is  excellent.  The  a^2  value  is  to  be  compared  with  measurements6 
near  the  zenith  ving  values  of  .  050  ±  .  019  nepers2  .  The  agreement  is 
hardly  outstanding,  but  seems  reasonable  considering  the  uncertainties  of 
the  quantities  involved  and  of  the  limited  size  of  the  samples  available. 
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Table  II.  1 

Bufton's  Flight  Data  --  Smoothed  and  averaged 

The  data  shown  here  is  based  on  smoothing  the  raw  data  over  1000  m 
altitude  increments.  This  is  accomplished  by  using  the  altitude  integrated 
measured  values  of  Q,2  .  We  plot  this  data,  passing  a  reasonably  smooth 

curve  through  the  results,  and  then  determine  the  change  in  the  integrated 
value,  as  represented  by  the  curve,  over  1000  m  increments.  Each  1000  m 
interval  is  centered  at  an  exact  number  of  km  altitude.  Thus  the  intervals 
are  0.  5  km  to  1.  5  km,  called  the  1  km  altitude  interval,  1.  5  km  to  2.  5  km, 
called  the  2  km  altitude  interval,  etc.  A  single  average  value  of  CN2  for 
each  interval  is  then  calculated  from  the  change  in  the  integrated  value. 


Alt.  (km)  Refractive-Index  Structure  Constant,  Q,2(h)x  1 016  ,  (m_21/3) 

_ I  Flt-  #4  j  Fit.  #6~|  Fit,  #7  |  Fit.  #9  !  Ave,  (4  Fits.  ) 

0.  47 


0.  15 


0.  19 
0.  23 
0.  17 


Fit.  #6 
0.  3  5 
1.  90 
0.  75 
0.  14 
0.  87 


0.  38 


0.  39 
1.  74 
0.  76 
1.  57 
0.  10 
0.  00 
0.  00 
0.  00 


0.  65 
1.  25 
0.  40 
0.  76 
0.  20 


0.  465 
1.  260 
0.  525 
0.  675 
0.  335 


Table  II.  2 

Button's  Flight  Data  --  Adjusted  Averages 


In  plotting  in  Fig.  U*  3  the  average  flight  data  from  Table  IL  1,  we  have 
noted  certain  fa'rly  obvious  anomalies.  We  believe  these  anomalies  are 
artifacts  due  to  the  limited  statistical  samplei  i.  e.  ,  the  use  of  aata  from 
only  four  flights.  These  anomalies  have  been  eliminated  in  the  adjusted 
averages  by  some  fairly  minor  redistribution  of  values  of  C^2  .  The 

adjusted  values  shown  in  this  table  correspond  to  the  curve  shown  in  Fig.  IL  3. 
An  asterisk  (*)  indicates  an  adjusted  value  that  deviates  from  the  average 
value  in  Table  II.  1. 


Altitude 
h  ,  (km) 

Refractive-Index  Structure  Constant 

CNS  X  1016  ,  (m-a*3) 

1 

0.  465 

2 

1.  260 

3 

0.  675  * 

4 

0.  525  * 

5 

0.  335 

6 

0.  190 

7 

0.  220 

8 

0.  250 

9 

0.  300 

10 

0.  395  * 

11 

0.  460  * 

12 

0.  510  * 

13 

0.  475  * 

14 

0.  385  * 

15 

0.  223 

16 

0.  153 

17 

0.  128 

18 

0.  090 

19 

0.  075 

20 

0.  080 
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Figure  II.  1.  Button's1  Thermal  Probe  Data,  Flight  #9 

The  values  of  CN2  computed  from  the  thermal  probe 
data  are  shown  without  any  altitude  smoothing.  The 
spikiness  of  the  data  is  so  severe  as  to  make  use  of 
the  data  for  calculation  of  atmospheric  optical  effects 
rather  difficult. 


Refractive-Index  Structure  Constant, 


i.1.0  I  XOI 


Figure  II.  3.  High  Altitude  Thermal  Probe  Measurement  of  the  Refractive -Index  Structure 
Constant.  The  data  is  extracted,  by  smoothing  and  adjusting  the  balloon  flight 
data  of  Bufton.1  The  circles  represent  the  averages  in  Table  H.  1.  and  the  curve 
represents  the  adjusted  data  in  Table  II.  2. 


I 


Theory, 


CHAPTER  III 

Differential  Angle  of  Arrival: 

Evaluation,  and  Measurement  Feasibility 


-  46  - 


Introduction 

In  this  chapter,  we  are  concerned  with  the  problem  of  developing 
theoretical  results  for  the  mean-square  difference  in  angle-of-arr ival  of 
light  seen  through  two  small  but  finite  apertures  viewing  the  same  point 
source.  The  two  apertures  are  separated  by  a  distance  which  is  a  few 
times  the  diameter  of  each  aperture.  We  wish  to  calculate  the  magnitude 
of  the  mean-square  angle  difference  as  a  function  of  the  separation. 

Our  interest  in  this  problem  is  related  to  an  experimental  problem. 
Ideally,  we  would  like  to  make  measurements  of  the  phase  difference  be¬ 
tween  the  two  apertures  when  viewing  a  monochromatic  point  source  (i.  e.  , 
a  laser).  This  would  give  us  a  direct  handle  on  measurement  of  the  phase 
structure  function,  which  because  of  its  relationship  to  the  atmospheric 
turbulence  limited  imaging  problem,  is  of  direct  interest  to  us.  For 
measurements  of  propagation  paths  between  two  points  on  the  ground, 
it  is  a  practical  matter  to  undertake  mean-square  phase  difference  meas¬ 
urements.  However,  when  the  source  is  off  the  ground,  as  it  must  be  if 
we  are  to  be  able  to  make  measurements  of  propagation  over  non-horizontal 
paths,  making  mean-square  phase  difference  measurements  is  not  practical. 
For  such  a  measurement,  the  receiver  will  have  to  track  the  source  and 
any  angular  tracking  error  (with  a  component  in  the  plane  defined  by  the 
line-of- sight  and  the  line-of- separation  of  the  two  apertures)  will  produce 
a  phase  difference.  This  phase  difference  will  appear  to  be  the  same  as 
a  propagation -induced  phase  difference.  Unless  the  tracking  is  virtually 
perfect,  the  measurement  results  will  be  dominated  by  this  tracking  error, 
leading  to  basically  spurious  conclusions. 

In  order  to  avoid  this  measurement  problem  and  yet  get  a  handle 
on  the  measurement  of  phase  differences,  one's  attention  is  then  naturally 
directed  to  the  possibilities  offered  by  making  angle -of-arrival  measure¬ 
ments.  The  angle  of  arrival  is  basically  the  first  derivative  of  the  phase 
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difference  and  we  should  be  able  to  extract  phase  difference  statistics 
information  from  angle -of-arrival  statistics. 

Actually,  the  basic  angle-of-arrival  measurements  are  plagued 
by  the  same  tracking  error  problem  as  is  the  phase-difference  measure¬ 
ments.  Tracking  error  effects  are  potentially  indistinguishable  from 
propagation-induced  angle-of-arrival  fluctuations.  However,  by  working 
with  the  difference  in  angle-of-arrival  at  two  separated  apertures  mounted 
on  the  same  tracking  assembly,  and  calculating  the  mean- square  value  of 
this  difference  (rather  than  the  covariance  of  angle-of-arrival  at  the  two 
apertures,  for  example),  we  avoid  any  effect  due  to  tracking  error.  The 
tracking  error  is  the  same  for  the  two  apertures,  and  its  contribution 
drops  out  of  the  difference  of  an'gle -of-arrival  at  the  two  apertures.  Hence 
our  interest  in  studying  the  expected  value  of  the  mean- square  difference 
of  angle-of-arrival  at  two  apertures  as  a  function  of  the  separation  of  the 
apertures.  This  problem  provides  the  theoretical  framework  for  a  prac¬ 
tical  measurement  of  effects  related  to  phase  distortion  over  a  non- 
horizontal  propagation  path. 

In  the  next  section,  we  define  the  measurement  problem  in  terms 
of  the  applicable  parameters  and  then  set  up  the  mathematical  formulation 
of  the  problem  in  a  tractable  form.  In  the  section  following  that,  v/e  carry 
out  the  necessary  analysis  to  reduce  our  formulation  to  a  two-dimensional 
definite  integral  which  can  be  evaluated  on  a  digital  computer.  In  the 
section  after  that,  we  present  the  computer  evaluation  procedure  and 
results.  The  final  section  presents  a  brief  discussion  of  these  results 
and  their  implication  in  terms  of  required  angular  measurement  precision. 

We  have  also  included  an  appendix  which  briefly  sets  forth  an 
alternate  approach  to  this  evaluation  problem  based  on  the  assumption 
that  the  apertures  are  very  small  and  that  the  angle-of-arrival  can  be 
considered  to  be  simply  the  derivative  of  the  phase.  We  show  that  this 
approach  runs  into  divergence  difficulties  which  can  be  lifted  by  introduction 
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of  the  inner  scale  of  turbulence.  We  argue,  however,  that  this  divergence 
is  more  properly  lifted  by  the  finite  size  of  the  measurement  aperture, 
and  that  a  theory  which  does  not  take  this  into  account  should  not  be  con¬ 
sidered  applicable  for  finite  sized  apertures. 

Problem  Definition  and  Formulation 

We  have  previously  shown1  that  over  a  circular  region  of  diameter 
D  ,  a  distorted  wavefront,  represented  by  the  random  function  0  (£)  ,  can 
be  conveniently  decomposed  into  a  set  of  terms  each  of  which  represents 
some  geometric  aspect  of  the  wavefront  distortion.  In  order  to  accomplish 
this  decomposition,  a  set  of  orthonormal  (two-dimensional  polynomial) 
functions,  Fj  (x)  ,  were  defined.  These  functions  were  directly  related 
to  the  Zernike  polynomials.  Of  particular  interest  to  us  here  are  the  two 
functions  which  are  related  to  tilt  along  the  two  orthogonal  component 
directions  in  the  x-plane.  If  we  denote  the  two  components  of  x  by  (x,  y) 
then  these  two  functions  are 


The  coefficient  a^^  is  obtained  from  the  equation 

^lit  =  J  dx  W(x,  D)  FtJIt  (x)  0(x)  ,  (2) 

where  the  integration  is  over  the  infinite  x-plane  with  the  actual  region 
of  integration  limited  by  the  circular  aperture  function  W(x,  D)  ,  defined 
by  the  equation 


W  (x,  D) 


1  if  |x|  s  £  D 

0  *  |x|  >  D 


(3) 


In  order  to  extract  the  tilt  angle,  a  from  atJlt  ,  we  need  a  scaling 
factor.  We  obtain  this  factor  by  noting  that  when  0(3c)  is  replaced  by  x 
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(or  y  ),  we  expect  the  slope  to  be  unity,  and  by  noting  that  the  slope 
associated  with  the  phase  function,  0  (x)  ,  should  be  calculated  from 
the  isophase  height  function  (\/ 2rr)  0  (x)  .  Since 

j d* w(*, D)  { (ibf x  }x  =  f  •  (4) 

it  follows  that  the  normalization  factor  is  (\/2rr)(nD4/64)  ^  .  This  means 
that  the  tilt  angle  a  associated  with  the  random  phase  function  0(x)  over 
a  circular  aperture  of  diameter  D  is 

a  =  Zn  (ttD4?  [  dx  W(x,  D)  Ft(lt  (x)  0(x)  .  (5) 

The  experimental  set-up  whose  results  we  wish  to  analyze  con¬ 
cerns  a  pair  of  circular  apertures  each  of  diameter  D  ,  with  their  centers 
at  Xj  and  3^  .  The  difference  of  the  angle- of- arrival  components  along 

the  x-axis,  as  seen  through  these  two  apertures,  can  be  written  in  accord¬ 
ance  with  Eq.  's  (1)  and  (5)  as 

O'!  "  Os  =  “^2  04"  J  dx  W(x,  D)  X  [0  (xj  +  x)  -  0  (Xj  +  x)]  .  (6) 

We  can  write  the  square  of  -  Cfe  as  the  product  of  two  integrals,  and 
that  as  a  double  integral  over  x  and  x'  .  Then  by  taking  advantage  of 
the  fact  that  we  can  commute  the  processes  of  integration  and  ensemble 
averaging  (which  we  denote  by  the  angle  brackets,  (  )  ),  we  see  that  we 
can  write  the  mean- square  angle- of-arr ival  x-component  difference  as 

<  (oa  -  Oa  )s  )  =  (-^i)8  JJ  dx  dx'  W £,  D)  W(x',  D)  x  x' 

X  <[0(xI  +  x)  -  <t>(xp  +  x)] [0(3^  +  xO  -  0(3fe  +  x')]>  .  (7) 

This  equation  represents  our  basic  problem  formulation.  In  the  next  sec¬ 
tion  we  shall  turn  our  attention  to  the  reduction  of  this  formulation  to  a 
numerically  evaluatable  result. 
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Formulation  Reduction 


At  this  point,  it  is  convenient  to  make  the  change  of  variables 
from  x  ,  x'  to  the  difference  and  sum  variables,  u  ,  v  ,  where 

U  =  X  -  x' 

V  =  £  (x  +  x') 

The  ensemble  average  in  Eq.  (7)  can  now  be  written  in  terms  of  the  four 
product  terms  expressed  as  phase  covariances,  as  (u)  -  (xg  - 3^  +  u) 

'  C0^‘*i  _u)  +  C0(u)  .  However,  by  adding  and  subtracting  twice  the 
phase  variance  and  appropriately  grouping  terms,  making  use  of  the  fact 
that  the  structure  function  for  a  stationary  random  variable  is  twice  the 
difference  of  variance  and  covariance,  we  can  more  conveniently  rewrite 
the  ensemble  average  in  Eq.  (7)  in  terms  of  this  phase  -  structure  function, 

as  +  *.00(Wu)  • 

If  we  introduce  the  vector  S  to  represent  the  center-tc- center 
separation  of  the  two  apertures,  i.  e.  , 


(8) 

(9) 


S  "  *3  -  *1 

then  we  can  rewrite  Eq.  (7)  as 

'32  X  \s 


(10) 


<fai-%)2)  =  //dSdv(v+ju),  (7-Ju),  W(7+J3,  D)  D) 

xt*^(S+u)+t^{S-u)-^(S)]  .  (11) 


In  the  above,  the  subscript  x  is  used  to  denote  that  only  the  component  along 
the  x-axis  is  to  be  considered. 


H  we  now  introduce  the  function  K(u,  D)  ,  where 

K(u,  D)  =  J  dv  tf+*u)x  W(v+fcu,  D)  W(7-*u,  D)  , 

then  we  can  rewrite  Eq.  (11)  as 


(12) 
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(<<H-aa»  =  (^r£)  J  du  K(3,D)[i^(S  +  -S)  +  i^(S  -  u)  -  i>0©]  .  (13) 

We  now  turn  our  attention  to  carrying  out  the  two-dimensional  integral 
in  Eq.  (12). 

To  carry  out  the  integral  in  Eq.  (12),  it  is  convenient  to  represent 
the  vector  v  by  the  two  components  (p,  q)  ,  where  the  p-axis  is  taken  to 
be  parallel  to  u  .  We  shall  let  6  denote  the  angle  between  the  x-axis 
and  u  .  We  can  then  write 


(v  ±ijf  u)x  =  p  cos  9  -  q  sin  0  ±  u  cos  0  .  (14) 

The  integration  limits  imposed  by  the  two  W-functions  in  Eq.  (12)  are  such 
that  we  can  now  write 

(  o- u)/  2  +  C<  V  z?-  (pW  8i3l1/Q 

K(u,  D)  -  2  J  dp  J  dq  £ps  cos2  q  _  2pq  cos  9  sin  0  +  q2  sin2  0 

-  (£  u)3  cos2  0  ]  (1  5) 

By  regrouping  and  adding  and  subtracting  various  quantities,  this  can  be 
rewritten  as 


Co-u)/2  +[(y  sf~  2)3]1/8 

K(u,  D)  =  2  j  dp  J  dq  [(p  +  fcu)3  cos2  0  -  (p  +  fcu)  u  cos3  9 

o  I* 

-C(e/3)a-(p-w/e)s] 

-  2  (p  +  £u)  q  cos  0  sin  9  +  q  u  cos  0  sin  9  +  q2  sin3  0]  (16) 

It  is  convenient  here  to  make  a  change  of  variable  so  that  (p  +£u)  is 
replaced  by  p  .  With  this  change  of  variable,  we  can  rewrite  Eq.  (16)  as 
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K(u,  D)  =  2  |*  dp  J  dq  (p2  cos2  Q  -  pu  cos3  9  -  2  pq  cos  9  sin  0 

U/3  -  0/3)3-pa]1/a 


+  qu  cos  0  sin  0  +  q2  sin  0) 


(17) 


Carrying  out  the  q- integration  is  a  trivial  proce s s.  We  obtain 

K(u,  D)  =  4  dp  [p3  [ftD)2-p2]*  cos3  0  -  p[ftD)2-p2]*  cos3  0 
a/3 

+  §  [  ft  D)2  -  P2  J3'*3  sin2  0}  .  (18) 

Now  if  we  extract  a  factor  of  ft  D)4  from  inside  the  integral  in 
Eq.  (18)  and  make  a  further  change  of  the  variable  of  integration,  replacing 
p/ ft  D)  by  p  ,  we  get 

K(u,  D)  =  4  ft  D)4  J*  dp  {p3  (1-p2)*  cos3  0  -  2p(l-p3)*(u/D)  cos3  0 

u/o 

+  5  (l-p2)3*5  sin3  0  }  .  (19) 

The  integrations  in  Eq.  (19)  can  be  carried  out  in  terms  of  the  formulas 
given  by  Dwight3.  We  can  write 

K(u.  D)  =  i  D4  ({i  cos'1  (u/D)  +  [1  -  (u/D)2]*[-i  (u/D)3  +  ±(u/D)]]  cos3  0 
-  2[  1  -  (u/  D)2  ]*[  -  §  (u/ D)2  +  -|]  (u/  D)  cos2  0 

+  {  8  c°s_1  (u/D)  +  [  1  -  (u/ Df  ]*  [_*_  (u/ D)3  - ~  (u/D)]}  sin3  0  ^  .  (20) 
This  can  be  rewritten  as 

K(u,  D)  =  A  D4 A  cos-1  (u/D)  +  [1-  (u/D)2]^-!-  (u/D)3  -_L  (u/D)] 

\  12  24 

+  [§  (u/D)3  -  |(u/D)]  [1  -  (u/D)2]*  cos2  0  }  )  (21) 
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If  we  now  combine  Eq.  's  (13)  and  (21),  we  obtain 


<  (cfc  -aa  )2  )  =  )4  ^  J  du  cos'1  (u/  D)  +  [  1  -  (u/  Df  ]*  {[£  (u/  Dp 

■  A  +  Es  (u/D)3-i  (u/D)]  cos2  9})[i^0  (s  +  u) 

+  *if0(S.^).^(S)]  .  (22) 


It  will  simplify  our  results  if  we  now  make  the  replacement 

J  =  S/D  ,  (23) 

and  after  bringing  in  a  factor  of  D-3  from  outside  the  integral,  make  a 
change  of  variables,  replacing  u/D  by  u  .  This  allows  us  to  rewrite 
Eq.  (22)  as 

<K-aa)2)  =  Jdu  cos-^u)  +  (1-u2)^  [(i  u3  -  JL  u) 

+  (5  u3  -  |u)  cos2  9]}{^0[D(^  +  u)] 

+  i*0[D(?  -  u)]  -J>0(Du)}  .  (24) 

At  this  point,  we  recall  that  we  define  9  as  the  angle  between  u 

and  the  x-axis.  Operationally,  the  x-axis  is  defined  by  the  fact  that  the 

component  of  angle-of- arrival  that  we  are  measuring  lies  in  the  x,  z-plane. 

We  will  pick  this  component  so  that  the  x-axis  makes  an  angle  \|t  with  the 

— * 

aperture  separation  vector  S  (and  with  3  )  .  This  means  that  the  angle 
between  u  and  J  is  0  +  ,  so  that  we  can  write 

|«/  ±  u|  =  [[of  ±  u  cos  (@  +  f)]2  +  [u  sin  (9  +  \|r)]a 

=  [/*  ±2^ucos  (9  +  t)  +  us]^  .  (25) 
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Taking  note  of  the  fact  that  the  pha se  -  structure  function,  ,  is  dependent 

on  only  the  magnitude  of  its  argument,  we  see  that  we  can  rewrite  Eq.  (24) 
as 

<  K  -  aa)2  )  =  (4-'\  A-J  j  de  Judu  (icos’Mu)  +  (1-u2)^  [(A  u3  -  4-  u) 

'  0  0  a4 

+  (|u3  -  ±u)  cos2  e]}  f  [D(y2  +  2^u  cos  (9  +  \|i)  +  u2]^} 

+  {  D[^2  -  2^u  cos  (0  +  t)+  u2  ]^)  -^0(°u)  )  •  (26) 

Eq.  (26)  is  the  most  general  expression  we  can  write  for  the  mean-square 
angle -of-arrival  difference.  Its  further  evaluation  can  be  carried  out  numer¬ 
ically  as  soon  as  we  have  an  expression  for  the  phase  structure  function,  . 

To  proceed  further  in  our  evaluation,  we  introduce  the  approximation 
that  the  phase  structure  function  can  be  accurately  represented  by  the  wave 
structure  function.  This  allows  us  to  write 

^0(x)  as  6.88  (x/r0)®/3  ,  (27) 

where  r0  is  a  length  representative  of  the  magnitude  of  the  wavefront 
distortion.  3  If  we  use  this  approximation  in  Eq.  (26),  we  obtain 

/4  .  -y  .1/3  .5/3  1 

<  (a,  -  a2)2  )  =  3.  44  (*)  (£)  (±  )  Jd0  Judu  {i  cos'Mu) 

v  0  o  o 

+  (1-u2)*  [(±  u3-  ^  u)  +  (i  u3-|  u)  cos2  0  ]} 

X  {[J2-  Zju  cos  (0+  f)+  u2]6/0+  [^2-2^u  cos  (0  +  f)  +  u2]®/®  -  2u®/3  } 

(28) 

In  this  form  our  expression  is  directly  suitable  for  computer  evaluation. 

We  take  up  the  matter  of  this  numerical  evaluation  and  presentation  of 
results  in  the  next  section. 
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Numerical  Evaluation 


To  proceed  beyond  Eq.  (28),  it  is  convenient  to  introduce  the  quan¬ 
tity  I^.ijr)  ,  which  we  define  by  the  equation 

2n 


I(0/,,\ji)  =  j  f  d9  J  du  u  jcos'1  (u) 


+  (1  -u2  u3  '  ~24  u)  +  ^  u3“  3  cosS  e  ]} 


X  ([  afi  +  2  JU  COS  (q  +  \J(  )  +  u2  ]B/6  +  [J2  -  Zj\l  cos  (9+ty )  +  u3  ]B/6  -  2  UB/3 

(29) 


This  allows  us  to  rewrite  Eq.  (28)  in  the  form 


<  (ax  -  a^)) 


3.  44 

Tt3 


(30) 


This  will  serve  as  our  final  result,  for  all  practical  purposes.  However, 
it  is  quite  interesting  to  take  note  of  the  expression  for  mean  square  angle 
of  arrival  variation  at  a  single  aperture,  (a2)  ,  as  given  in  Appendix  B, 

Eq.  (B-2),  and  recognize  that  we  could  rewrite  Eq.  (30)  as 

<(oa  -%?)  ~  <a3>  I U,  Mr)  .  (31) 

Since  o  is  a  two-component  tilt,  i.  e.  ,  it  includes  both  x  and  y  components 
in  the  mean  square,  while  -  Og  represents  the  difference  of  two  single 
component  tilts,  it  would  be  necessary  for  Oj  and  a2  to  be  anti-correlated 
for  1{J,  t)  to  have  a  value  greater  than  unity,  and  in  no  case  could  it  be 
greater  than  two.  In  fact,  we  find  it  never  gets  to  be  quite  as  large  as  unity. 
The  numerical  evaluation  of  It/,i)f)  is  a  straightforward  computer  program¬ 
ming  problem.  We  have  prepared  a  program  to  carry  out  these  calculations. 
The  program  is  listed  in  Table  III.  1,  and  the  results  for  I(^,  0)  and  H^.^tt) 
are  given  in  Table  III.  2  for  S=  1  to  20  in  steps  of  ^  .  As  can  be  seen  from 


these  results,  for  =  0  ,  I  has  values  in  ihe  range  of  0.  5  to  0.75, 
approximately,  while  for  t|f  =  ^  tt  ,  the  values  are  somewhat  smaller. 

Discussion  of  Results 

Our  basic  interest  in  all  of  this  work  has  been  in  assessing  the 
practicality  of  using  mean  square  difference  in  angle  of  arrival  measure¬ 
ments  to  determine  the  optical  strength  of  turbulence  along  the  propagation 
path,  or  more  succintly  put,  to  determine  rQ  .  The  critical  question  is 
how  much  larger  <  (o^  -  Og)2  )  is  expected  to  be  than  our  measurement 
precision.  From  Eq.  (30)  and  the  fact  that  Hj,  ill)  has  a  value  of  the 
order  of  0.  5  ,  we  see  that  the  mean  square  measurement  value,  which 
we  call  the  signal  squared,  should  be 


(Sig)2 


3.  44  ,X'1/3  /_Lf/3 

ns  \rQ  ) 


(32) 


Our  equipment  will  make  measurements  of  single  ax's  angle  of  arrival  with 
some  rms  uncertainty,  69  .  The  mean  square  variability  in  measure¬ 
ment  of  (  (otj  -a2)2)  1  which  is  our  noise  squared,  is  just  the  sum  of  the 

error  for  and  for  a B  .  Thus  we  can  write  for  the  noise, 


(Noise)2  as  2  (69)2 


(33) 


The  rms  error  in  determination  of  rQ 
be  determined  from  the  equation 

3  44 

(Sig  +  Noise)3  «  “ 

Solving  for  6r0  from  Eq.  's  (32), 
we  get 


,  which  we  denote  by  6r0  ,  can 


(33),  and  (34),  using  6r0  «  r0  , 
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2.87  oe 


(35) 


We  see,  for  example,  that  if  we  work  at  X  =  0.  633  p,m  ,  use  aperture 
diameters  of  D  =  0.  1  m  ,  and  if  rQ  m  0.  1  m  ,  (for  which  rms  angle  of 
arrival  variation  is  (3.44/n 2)1^2  (X /  D)1^6  (X/tq)6^6  =  3.  74  p,rad  ),  then  we 
can  achieve  a  measurement  precision  6  r0  / r0  «  4.  53  X  106  59  .  To  get 
a  10%  accuracy,  the  rms  ei.’or  60  would  have  to  be  no  greater  than  0.22 
p,rad  =  0.  045  arc  seconds.  If  69  »  0.  1  arc  seconds,  then  6r0  /rQ  would 
be  of  the  order  of  22%.  If  the  propagation  path  were  such  that  r0  were 
0.  5  m  (for  which  the  rms  angle  of  arrival  variation  is  6.  66  p,rad),  then 
we  would  have  6r0  / r0  »  2.  54  X  1  6  9  •  To  get  a  10%  accuracy,  the 

rms  error  in  69  would  have  to  be  no  greater  than  0.  39  u,rad  =  0.  081  arc 

seconds. 

The  feasibility  of  the  experiment  obviously  depends  critically  on 
our  ability  to  select  propagation  paths  for  which  r0  will  be  reasonably 
small,  and  on  our  ability  to  achieve  a  small  enough  value  of  60  . 

If  we  can  reliably  obtain  69  «  0.  1  arc  seconds,  it  should  be  possible 

to  obtain  viseful  measurement  results.  The  basis  for  the  data  reduction 
of  any  such  results  is  provided  by  Eq.  (30),  and  the  data  in  Table  III.  2. 
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APPENDIX  A  of  Chapter  III 

An  alternate  approach  to  the  problem  of  treating  angle-of-arrival 
statistics  is  to  use  the  derivative  of  the  phase  function  as  a  measure  of 
the  angle  of  arrival.  Quite  obviously,  the  derivative  represents  the  wave- 
front  tilt  in  a  region  of  zero  extent,  and  at  least  nominally  there  should 
be  no  difficulty  in  utilizing  such  an  approach.  We  shall  see,  however, 
that  this  approach  involves  us  in  a  divergence  problem.  With  this  formu¬ 
lation,  the  divergence  is  lifted  by  recourse  to  the  inner  scale  of  turbulence, 
providing  smoothing  over  what  otherwise  would  be  a  cusp  in  the  statistics, 
i.  e.  ,  x^3  has  i  cusp  at  x  =  0  .  In  an  experiment,  the  smoothing  would 

be  provided  by  the  finite  size  of  the  aperture  diameter  D  ,  and  results 
obtained  using  the  derivative  approach  will  depend  on  the  inner  scale  of 
turbulence  JlQ  ,  where  they  ought  to  show  dependence  on  D  . 

Using  the  derivative,  we  can  write  for  the  x- component  of  the  local 
angle-of-arrival  at  x  , 


so  that  the  difference  in  x-components  of  the  local  angle-of-arrival  at 
Xj  and  Xg  can  be  written  as 


ax  ‘  °a 


#(S»>  -  *<*>] 


(A-2) 


The  mean  square  difference  can  then  be  written  as 

«5>  '  *W>  ■  (•£)’  <  [£#*>  - 


*3  =*8 


(A-3) 


(This  rather  unusual  notation  introducing  x^  and  x^  is  necessary  so  that 
when  we  carry  out  the  multiplication  and  rearrange  the  terms,  we  can  tell 
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which  differentiation  goes  with  which  phase  function.  )  Carrying  out  the 
multiplication,  rearranging  terms,  and  interchanging  the  order  rf  ensemble 
averaging  and  differentiation,  we  obtain 


U^i-Sa)3)  =  (2Xn)  ^  <0(^)0(3^)>  ^  <0(^1  0(%» 


^  <0  ^  0  +  aL,  <0  ^  0  Xj  (A-4) 


9  a 


Xx  =  Xi 
-4*  ->* 

xa  =  xa 


The  ensemble  averages  in  Eq.  (A-4)  can  be  recognized  as  the  phase 
covariance  function.  By  adding-and- subtracting  the  phase  variance,  we  can 
obtain  the  phase  structure  function,  so  that  Eq.  (A-4)  can  be  rewritten  as 


«5> -  W2>  ■  (£)’  [4 £ 4 ** -*>  +  4^  £ 


+  ^i^,  ' 4  8Xj  ^  •  (A_5) 

xa  =  xa 

We  now  note  that  if  we  introduce  appropriate  sum  and  difference  coordinates, 
and  hold  the  sum  constant,  then 


9x,  9Xj  3(t>  "  X1 ) 


^ 


(A-6) 


where  u  =  xt  -  ,  and  ux  denotes  the  x-component  of  u  .  This  means 

that  we  can  rewrite  Eq.  (A- 5)  as 


<<s,-5,)s>=-(^)2  {^-ys>  -  [^  •V1')]3.=0}  •  <A-7> 


where 


S  = 


(A-8) 
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If  we  substitute  the  approximate  form  for  the  phase  structure 
from  Eq.  (27)  of  the  main  text,  the  nature  of  our  divergence  problem 
becomes  manifest.  We  obtain 

<6,  -S,>2>=-(^)2  6.  88  r0-e/3{l  <S„=+  S,*)'™  (§  S,=  +  S,=) 

-  f|-(s.'2+  s;2)-"a(§sT'a+  s/s>]s-=0  }  •  (A'9) 

s;=o 

— ♦ 

The  second  group  of  terms  in  Eq.  (A-9)  is  clearly  divergent  when  S'  is 
set  equal  to  zero.  We  believe  that  the  first  group  of  terms  accurately 
represents  the  S-dependence  of  the  mean  square  difference  of  the 
x- component  of  angle -of- arrival,  but  to  get  any  use  out  of  the  result,  we 
must  find  a  way  of  renormalizing  the  divergence  in  Eq.  (A-9).  This  re¬ 
normalization  can  be  provided  by  considering  that  Eq.  (27)  is  only  valid 
for  values  of  x  greater  than  or  equal  to  the  inner  scale  of  turbulence  Iq  . 
We  argue  that  more  properly  the  divergent  term  in  Eq.  (A-9)  should  be 
represented  by  the  derivative  evaluated  at  l0  rather  than  at  zero.  This 
then  leads  to  a  finite  term  of  the  order  of  £o~1/3  in  place  of  the  divergent 
term,  and  a  result  similar  to  our  exact  one  obtained  in  the  main  text,  but 
with  the  role  of  the  aperture  diameter  D  filled  by  the  inner  scale  of  tur¬ 
bulence,  SLq  .  Since  for  any  practical  experiment  we  will  make  D  greater 
than  lQ  ,  we  conclude  that  the  applicable  renormalization  would  involve 
D  rather  than  A0  ,  with  the  accurate  result  being  that  obtained  in  the 
main  text  using  Eq.  (2)  to  obtain  the  tilt  rather  than  the  derivative  pro¬ 
cedure  of  Eq.  (A-l)  in  this  appendix. 
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APPENDIX  B  of  Chapter  III 

It  is  necessary  in  the  text  to  reference  our  mean  square  difference 
in  angle  -  of -arrival  results  to  the  mean  square  angular  variation  for  meas¬ 
urement  on  a  single  aperture.  This  calculation  of  the  mean  square  angular 
variat  on  has  been  carried  out  in  reference  1,  but  without  normalization. 

In  this  appendix,  we  shall  use  those  results  and  only  concern  ourselves 
with  the  normalization. 

If  we  let  ax  and  ay  denote  the  x  and  y  components  of  tilt  cor¬ 
responding  to  amt  defined  in  Eq.  (2),  then  from  Eq.  (7.8a)  of  reference  1*. 
we  have 


((a,)3  +  (ay  )2 )  =  0.883  (In  CP)  (D/r^  .  (B-l) 

As  pointed  out  in  the  discussion  just  after  Eq.  (4)  in  the  text  of  this  report, 
a  factor  of  (X / 2rr )(rr  D4/64)'VE  is  required  to  convert  the  atm  coefficient 
to  an  angle-of-arrival  coefficient  a  .  Thus  the  mean  square  angle-of-arrival 
fluctuation  for  a  single  aperture  of  diameter  D  ,  summing  the  two  com¬ 
ponents,  should  be 


<a2>  =  <(ax)2  +  (a^)2) 


*  There  is  a  factor  of  in  D2  missing  from  the  right-hand-side  of  Eq.  's 
(7.  8a,  b,  and  c)  of  reference  1.  This  can  be  traced  back  to  a  failure  to 
pick  up  a  factor  of  ttrs  in  solving  Eq.  (4.  4  ')  to  obtain  Eq.  's  (4.  6a,  b, 
and  c)  in  that  paper. 
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Table  III.  2 
Calculated  Values 
£°r  !(-/  ,  ♦  ) 


J 

1(^.0) 

*  (^>£17) 

1 .  P 

, +4. 599E-P 1 

+9. 2E4E-P 1 

1  .  5 

+5. 4537-P 1 
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9.7 

+5. 9P7E-P 1 

+3. 0 1 5E-P 1 
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